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O ' 1. Introduction 



ABSTRACT 



Aims. The radiation from newly born stars photoevaporates their parental neutral cloud, leading to the formation of dense clumps that 
will eventually form stars. 

Methods. We present 3D simulations of the interaction of a neutral cloud with an external ionising radiation field, and compute the 
angular momenta of these collapsing clumps. 

Results. The angular momenta of these collapsing clumps show that they have preferential orient mostly perpendicular to the direction 
of the incident ionising photon field. Therefore, the axes of the jet systems that will be eventually ejected (from the star + accretion 
disk systems that will form) will be oriented approximately perpendicular to the direction to the photoionising source. 

Key words. ISM: kinematics and dynamics - ISM: clouds - ISM: HII regions - stars: formation - ISM: individual objects: HH 555, 
HH 666, HH 333 
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The radiation from newly born stars photoionises and erodes the 
parental cloud, producing structures such as the so-called ele- 
phant trunks. At the head of an elephant trunk, the interaction 
of the shock (driven by the photoevaporation process) with pre- 
viously existing density perturbations leads to the formation of 
dense clumps. Some of these clumps might have enough mass 
to be self-gravitating, and will eventually form young stars that 
eject bipolar outflows. We describe observed examples of this 
kind of configuration. 

Bally & Reipurth (2003) discovered HH objects in the 
molecular cloud associated with the Pelican Nebula, including 
HH 555. This outflow emerges from the tip of a long elephant 
trunk, providing direct evidence of ongoing star formation in this 
region. The outflow axis of HH 555 is approximately perpendic- 
ular to the elephant trunk, which is aligned with the direction to 
the photoionising source. 

Another example of this kind of configuration is the HH 666 
outflow in the Carina nebula. HH 666 also emerges from close 
to the tip of an elephant trunk, and its axis is almost perpendic- 
ular to the direction towards ri Carinae (Smith, Bally & Brooks 

2004) . An HST image of this region (Bally, Reipurth & Davis 
2007), shows a second jet emerging from a nearby elephant 
trunk, with a direction almost parallel to the HH 666 outflow. 

A final example is provided by HH 333. This jet emerges 
from the tip of an elephant trunk within the complex Trifid neb- 
ula (M20). It is a single- sided jet with measured radial veloci- 
ties (Rosado et al. 1999) and proper motions (Yusef-Zadeh et al. 

2005) that indicate it has the kinematical properties of a standard 
HH jet. Again HH 333 has an outflow direction approximately 
perpendicular to the direction to the ionising source. 

Reach et al. (2009) presented observations of the elephant 
trunk of globule IC 1396A. They detected outflow activity from 
a number of young stars in the region. However, it is impossible 
to determine the outflow axes from these observations. 



Even though the number of four outflows (two in the HH 666 
region, see above) observed to be emerging from tips of elephant 
trunks is quite small, the alignment approximately perpendicular 
to the direction to the ionising photon source might be indica- 
tive of a systematic alignment. This alignment implies that the 
angular momenta of the low mass star-Fdisk systems producing 
outflows from stellar sources in the tip of elephant trunks are 
more or less perpendicular to the direction of the ionising pho- 
ton field (produced by the massive stars giving rise to the pho- 
toionised nebulae and elephant trunk structures). These angular 
momenta presumably preserve the direction of the rotation axes 
of the dense clumps that collapsed to form the outflow sources. 

In the present paper, we explore the interaction between an 
ionising photon field and an environment with density perturba- 
tions. This interaction produces elongated structures reminiscent 
of elephant trunks, with dense, embedded clumps. In particular, 
we focus on whether or not these dense clumps have angular 
momenta preferentially oriented perpendicular to the direction 
towards the photoionising source. 

Mellema et al. (2006) carried out 3D, radiation gasdynamic 
simulations of an H II region expanding in an ISM with power- 
law density perturbations. They find that this configuration nat- 
urally leads to the formation of dense, radially elongated struc- 
tures, which resemble elephant trunks. Also, Gahm et al. (2006) 
study the role of magnetic fields in the formation of elephant 
trunks. Finally, Gritschneder et al. (2009) carried out a simula- 
tion of an initially plane ionising front travelling into a structured 
medium. 

Our work emulates the approach of Mellema et al. (2006) 
and Gritschneder et al. (2009). We focus on a small region of the 
edge of an expanding H II region, and carry out a 3D radiation 
gasdynamic simulation (including the self-gravity of the gas) of 
the formation of a dense, neutral structure. We then identify high 
density clumps within this "elephant trunk", and compute their 
angular momenta. Finally, we study the mass distribution of the 
clumps, and the distributions of the orientation and magnitude 
of their angular momenta. 
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The paper is organized as follows. In Sect. 2, we describe 
the gasdynamic code and the parameters used for the numerical 
simulation. The results from the simulation and the clump statis- 
tics are presented in Sect. 3. Finally, our results are summarised 
in Sect. 4. 

2. Code & settings 

2.1. Code 

We carried out a 3D simulation with a code that solves the 3D 
gasdynamic equations, the Poisson equation for the gravitational 
field, and a rate equation for neutral hydrogen, including the 
transfer of ionising photons at the Lyman limit. The gas is ini- 
tially atomic, and the models do not consider the photodissoci- 
ation of molecular material because of the presence of a FUV 
radiation field. This code was described by Raga et al. (2008). 

We modified the code of Raga et al. (2008) to include the 
"two temperature" equation of state described by Esquivel & 
Raga (2007, hereafter E07). This equation of state assigns tem- 
peratures between 10 K (for neutral gas) and 10^ K (for gas 
with fully ionised H) with a linear dependence on the H ioni- 
sation fraction. Therefore, instead of solving an energy equation 
with the appropriate heating and cooling terms (see Raga et al. 
2008), we replace it with this two-temperature equation of state. 
We also included the self-gravity of the gas. We use a succes- 
sive over relaxation (SOR) method to solve the Poisson equation 
for the gravitational potential, and then include the gravitational 
force in the momentum and energy equations. We do not include 
a treatment of the diff'use, ionising photon field. 

2.2. Settings 

The computational domain has a size of (3.0, 1.5, 1.5) x 10^^ cm 
(along the x-, y- and z-axes, respectively), which is resolved with 
a uniform grid of 256x 128x 128 grid points. We impose transmi- 
sion boundaries in the x-direction and periodic boundaries along 
the y and z-directions. The periodic conditions are imposed in 
the gasdynamic equations, in Poisson's equation (for the gravi- 
tational field), and in the radiative transfer equations. 

We start with an inhomogeneous density structure with a 
power-law power- spectrum index of -11/3 (i.e. P[k] oc k~^^^^, 
where k is the wave-number), as described in Esquivel et al. 
(2003). The initial density structure does not have any motion. 
To simulate the edge of an HII region, the computational do- 
main is divided into two portions with a dividing line at x = 
4 X 10^^ cm from the left edge of the domain. The portion to 
the left is filled with an ionised medium (with a temperature of 
10"^ K, and the portion to the right is filled with a neutral medium 
(with a temperature of 10 K). The average density in the neutral 
medium is a factor of 100 higher than the one in the ionised 
medium, and the transition between the two (also in terms of 
temperature and ionisation fraction) follows a tanh profile with a 
width of ~ 10 pixels. The resulting neutral structure has a mass 
of 228 Mq. 

To calculate the gravitational field, we only consider the 
gravitational force resulting from the density perturbations. In 
other words, we subtract a density po = 3.5 1 x 10"^^ g cm"^ (cor- 
responding to the lower density regions in the initial distribution 
of neutral material) from the density used in Poisson's equation. 
In this way, we avoid a generalized collapse of the dense slab 
structure that fills the computational domain. We also run a sim- 
ulation in which the gravitational force was "turned off'" to illus- 
trate the eff'ect of the self-gravity of the gas. 
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Fig. 1. Time evolution of the xy mid-plane density stratification 
without self-gravity. The three frames are labeled with the cor- 
responding elapsed times. The density stratifications are shown 
with the logarithmic grey scale given (in g cm"^) by the top right 
bar. In the three frames we also show the contour corresponding 
to an H ionisation fraction of 50%, which indicates the position 
of the ionisation front. The x and j-axes are labeled in cm. 



A plane-parallel ionising photon field Fq = 8.8 x lO^^cm^s"^ 
is incident on the computational domain along the x-axis. This 
photon ffux corresponds to a star with an ionising photon rate 

= 10"^^ s"^ located at a distance D = 9.5x 10^^ cm from the 
edge of the computational domain in the -x direction. 



3. Results 

We allowed the model to run from the initial conditions de- 
scribed in Sect. 2, to a ^ = 200 kyr evolutionary time. Figure 
1 shows the time evolution in the mid-plane density stratifica- 
tion without including the self gravity of the gas. Figure 2 shows 
the same simulation but adding the force that arises from self- 
gravity. 

From both Figs., it is evident that the ionisation front be- 
comes highly corrugated with dense condensations at the tip of a 
number of protruding "fingers". At ^ = 100 kyr, the eff'ect of self- 
gravity is only to produce denser condensations at the tip of the 
fingers. At ^ = 200 kyr, however, the density structures obtained 
without (Fig. 1) and with self-gravity (Fig. 2) are quite diff'erent. 
In the self-gravitating simulation, a dense, central structure (de- 
tached from the ionisation front and absent in the non-gravitating 
simulation) is produced. 

Given the important diff'erences found when including self- 
gravity, we present an analysis of clump formation only for the 
ffow obtained from the self-gravitating simulation. Interestingly, 
if one repeats the analysis for the non-gravitating simulation, 
similar results are found (these results are not shown in the 
present paper). 
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Fig. 2. Same as Fig. 1, but for the simulation that includes the 
self-gravity of the gas. 



- p,= 10-20 




50 100 150 

time [kyr] 



Fig. 3. Number of neutral clumps as a function of time, obtained 
for three different density cuttoffs. The results correspond to the 
simulation which includes self-gravity (see Fig. 2). 



To quantify the number of clumps produced, following E07 
we calculate the number of spatially connected neutral struc- 
tures with densities above a specified cutoff' density pc. In par- 
ticular, we choose cutoff' values of pc = 3x, 10"^^, 10"^^ and 
10-2^ gcm-^ 

The number of clumps (fragments) obtained for diff'erent 
density cutoff's, is shown as a function of time in Fig. 3. To 
determine the number of clumps, we consider time intervals of 
20 kyr (corresponding to the width of the bins in the histograms 
of Fig. 3). We then calculate the number of clumps in 8 outputs 
within each of these time intervals, and then compute an average 
number of clumps for each time interval. 

For the lowest cutoff' density (pc = 10"^^), the initial density 
distribution has ~ 80 clumps, and the number of clumps first de- 
creases with time, stabilizes at ~ 5 for 40 < ^ < 170 kyr, and 
then continues to increase a little ait > 170 kyr (see Fig. 3). For 
the intermediate cutoff' density (pc = 10"^^), the initial distribu- 
tion has ~ 500 clumps and the number of clumps first decreases 




Fig. 4. Fraction of the clumps (obtained for a cutoff' density pc = 
3x 10"^^g cm"^) with diff'erent orientations a (between the x-axis 
and the xj-projections of the angular momenta of the clumps). 
The panels are labeled with the elapsed time corresponding to 
the time-frames from which the three angular distributions were 
obtained. The results correspond to the simulation that includes 
self-gravity (see Fig. 2). 

and then remains approximately constant as a function of time 
(with a value of ~ 20). 

The initial density distribution has no clumps with densities 
above the highest chosen cutoff' density, pc = 3 x 10"^^ g cm"^ 
(see above). Interconnected structures of sufficiently density 
only start to appear at ^ ^ 30 kyr, and their number grows 
monotonically with time, stabilizing at a number of ~ 20 for 
^ > 110 kyr (see Fig. 3). For each of the detected clumps, we 
first compute the position of the centre of mass : 



(1) 



M Jy 

where V is the contiguous volume of the clump and 

M= r pJ^x, (2) 

is its mass. We then compute the angular momentum with re- 
spect to the centre of mass of each clump 



L = J"p (r - Rcm) X V J^x . 



(3) 



We assume that we observe the computed ffow along the z-axis 
(i. e., that the xj-plane of the computational domain is parallel 
to the plane of the sky). The angle 



a = I SLTCtan Ly / Lx\ 



(4) 



(with L given by equation [3]) then corresponds to the orientation 
angle (with respect to the direction of the ionising photon field) 
of the angular momentum of the clumps projected onto the plane 
of the sky. As discussed in Sect. 1 , these directions correspond 
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Fig. 5. The top panel shows the specific angular momentum 
(L/M) of each clump as a function of the clump's mass. The 
black line shows the angular momentum associated with the 
outer orbit of an accretion disc radius of = 50 AU around 
a M = 2Mq star (L/M ~ 5 x 10^^ km^s'^ see the text). The 
bottom panel shows the orientation a (between the x-axis and 
the xy-projection of the angular momentum) for each clump, as 
a function of its mass. The triangles show the clumps found for 
an elapsed time of 95 kyr and the crosses for an elapsed time of 
195 kyr. The clumps were obtained using apc = 3 x 10"^^g cm"^ 
density cutoff, using the results from the simulation of Fig. 2. 



to the directions in which bipolar outflows will eventually be 
ejected when (and if) the clumps form star + accretion disk sys- 
tems. 

In Fig. 4, we show histograms indicating the fraction of 
clumps (obtained for a cutoff density = 3 x 10~^^g cm~^) with 
different orientations or, for the three different elapsed times (t = 
95, 145, and 200 kyr). For early times, we find that the a values 
of the clumps are randomly distributed (between ~ 40 and 180°). 
For t = 200 kyr, « 36% of the clumps have 70° < or < 100°, and 
more than ^ 55% of the clumps have 60° < a < 100°. From this 
result, we conclude that the dense clumps being formed have 
angular momenta preferentially aligned in directions perpendic- 
ular to the direction of the incident ionising photon field (which 
is parallel to the x-axis). 

The bottom panel of Fig. 5 shows the projected orientation 
a of the angular momentum as a function of clump mass for all 
of the clumps obtained with the pc = 3 x 10" ^^g cm"^ cutoff 
density, for elapsed times t = 95 and 200 kyr. We see that at 
t = 95 kyr most of the clumps (triangles) have masses Mc < 
0.3Mo and angular momenta with all a orientations. For t = 
200 kyr (crosses, see bottom panel of Fig. 5), we see that all of 
the clumps with IM© < Mc < 60Mo have angular momenta 
with orientation angles 70° < a < 120°. The lower mass clumps 
(with Mc < 0.3Mo) have angular momenta with more widely 
distributed orientations. 

We computed the moduli of the specific angular momenta 
(i. e., the momentum per unit mass) of the clumps. The values 
of L/M for t = 95 and 195 kyr are shown as a function of the 
clump mass M in the top panel of Fig. 5. For ^ = 95 kyr, we see 
that the clumps with Mc < 03Mq have M/L < 6 x 10^^ km^s"^. 
The more massive clumps (with Mc > IM© ) have M/L > 6 x 
10^^ km^s-^ For t = 200 kyr, we see that the Mc < O.O2M0 
clumps have M/L < 6 x 10^^ km^s"\ while the clumps with 
Mc > O.O2M0 have M/L > 6 x 10^^ km^s"^. 
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Fig. 6. Time evolution in the most massive, neutral clump ob- 
tained from the simulation with self-gravity (see Fig. 2). The top 
panel shows the mass, and the central panel shows the ratio of 
the clump mass to the Jeans mass as a function of time. The bot- 
tom panel shows the time evolution in orientation a of the angu- 
lar momentum of this clump (where a is the angle between the 
X-axis and the x}^-plane projection of the angular momentum). 



We now evaluate whether or not these specific angular mo- 
menta have values comparable to those observed in young star 
systems. Typical T Tauri stars have masses M » 2Mq and ac- 
cretion disks with radii rj) ~ 50 AU. The outer Keplerian orbit 
of the disk then has a specific angular momentum (L/M)d ~ 
5 X 10^^ km^s"^. This outer orbit is determined by the material 
of the highest angular momentum in the core from which the 
star-Fdisk system was formed (see, e. g., Ulrich 1976). 

This value of (L/M)d is shown with a horizontal line in the 
top panel of Fig. 5. It is clear that many of the clumps formed 
in our simulation have specific angular momenta that are sub- 
stantially higher than deduced from the radius of disks around 
T Tauri stars. From this, we conclude that the angular momenta 
of the clumps generated in our simulation are substantial. 

A relevant question is whether the clumps obtained in our 
simulations are resolved well enough for the calculations of an- 
gular momenta to be meaningful. As an example, we consider 
the clumps found for the pc = 3 x 10"^^ g cm"^ cutoff density 
at time t = 200 kyr. The lower mass clumps (see Fig. 5), of 
M ^ 0.01 Mo, are resolved with ~ 5 grid points. The clumps 
of M ^ 0.1 M© are resolved with ~ 50 grid points. The clumps 
with 1 < M < 100 Mo are resolved with ~ 500 to 5000 grid 
points. Therefore, for clumps with M > 0.1 Mo, the resolution 
of the internal structure of the clumps (with 50 grid points, cor- 
responding to ~ 4 grid points along each axis) appears to be 
appropriate for obtaining a meaningful estimate of the angular 
momentum. 

From the number of grid points Nc within the clumps, we 
can estimate the characteristic radii of the clumps to be Rc = 
0.012 pc X A^y^/2, where 0.012 pc is the resolution of the com- 
putational cells in our numerical simulation. From the values of 
Nc given in the previous paragraph, we then see that the clumps 
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Fig. 7. Density stratification and flow field in a region around 
the centre of mass of the most massive clump at t = 200 kyr 
(see Fig. 6). The top panel shows the flow on a xy-cut and the 
bottom panel the flow on a xz-cut The origin of the coordinate 
system coincides with the centre of mass of the clump (obtained 
with a pc = 3 X 10"^^ g cm"^ density cutofl', which corresponds 
to a number density of 1.4 x 10^ cm"^). The colour scale shows 
the logarithmic density distribution (given in cm"^ by the bar 
on the top right). The arrows show the flow velocities on the 
planes of the two cuts for densities p > 10"^^ g cm"^. An arrow 
with a length equal to the separation between successive arrows 
corresponds to a flow velocity of 0.5 km s"^ . 

obtained from our simulation have characteristic radii Rc ~ 0.01, 
0.02, 0.05, and 0.1 pc, for clump masses of M ^ 0.01, 0.1, 1, and 
100 Mq, respectively. 

Finally, we study the evolution in the most massive, neutral 
clump (detected with p^ = 3 x 10"^^ g cm"^). As seen in Fig. 6, 
this clump has a mass that grows monotonically from 0.6 M© at 
^ = 38 kyr, to 60 M© at ^ = 200 kyr. The orientation angle a (on 
the plane of the sky) of its angular momentum stabilizes rapidly 
at a ^ 90° for t > 70 kyr. 

We compute the Jeans mass of this clump to 



(see, e. g., E07) where G is Newton's constant, Cs is the sound 
speed of the neutral medium (see Sect. 2), and p = MdV^^^. We 
show the ratio MdMj for the most massive clump as a function 
of time in the central panel of Fig. 6. It is clear that this clump is 
Jeans unstable for -160 kyr, which is a long enough timescale 
for the formation of a low mass star. 

Figure 7 shows the density and flow velocity distributions in 
the xy- and xz-planes, within a (2 x 10^^ cm)^ region centred on 
the centre of mass of the most massive clump at ^ = 200 kyr. In 
the two cuts that are shown, we see that the region with densities 



higher than the pc = 3 x 10"^^ g cm"^ density cutofl' (which 
corresponds to a number density of 1 .4 x 10^ cm"^) has a number 
of density maxima, none of which coincides with the centre of 
mass of the structure. The xz-plane (bottom panel) shows the 
velocity field that gives rise to the angular momentum of the 
clump. 

4. Conclusions 

We have presented the results of numerical simulations of a neu- 
tral structure with power-law density perturbations that is photo- 
evaporated by an incident, plane-parallel, ionising photon field, 
with and without the self gravity of the gas. In this interaction, a 
number of dense, neutral clumps are produced. Our simulations 
are similar to those presented by Gritschneder et al. (2009). The 
main diff'erence is that while they started their simulations in a 
medium with turbulent motions, our simulations begin in a sta- 
tionary medium with density perturbations. In our simulations, 
the velocity field that develops is therefore mainly the result of 
the interaction with the ionising photon field. Defining clumps 
as contiguous structures above a cutoff' density pc, we compute 
the statistics of the number of clumps as a function of elapsed 
time (for diff'erent values of pc). We then fix the cutoff' density at 
Pc = 3 X 10"^^ g cm"^, to focus on the denser clumps appearing 
at later elapsed times. 

For these clumps, we compute the vector angular momenta, 
from which we obtain the direction of the rotation axes (pro- 
jected on the plane of the sky) and the specific angular momenta. 
We find that as a function of evolutionary time we obtain orienta- 
tions that are aligned increasingly perpendicular to the direction 
of the incident, ionising photon field. 

For the most massive clump, we find that it has a mass that 
increases across the range ^ 0.6-60 M© (during the ^ = 38 ^ 
200 kyr period), and that the orientation angle a of its angular 
momentum eventually stabilizes at a ^ 90° (i. e., the direction 
perpendicular to the direction of the incident photon field). We 
use an estimate of the Jeans mass of the clump to show that it 
is Jeans unstable throughout the ^ = 38 ^ 200 kyr period. This 
timespan is long enough for a low mass star to form within the 
most massive clump. However, at the resolution of our simu- 
lation (with a grid spacing of ^ 800 AU), we naturally do not 
succeed in form a star-hdisk system. 

If we analyse our non-gravitating simulation (see Fig. 1), we 
obtain qualitatively similar results. Regardless of whether we 
consider the self-gravity of the gas or not, we produce clumps 
with angular momenta preferentially aligned perpendicular to 
the direction of the incident ionising photon field. Even though 
it is impossible to provide a full explanation of this alignment, 
it is possible to provide a qualitative explanation. During the in- 
teraction of the ionising photon field with a perturbed density 
structure, a corrugated ionisation front is produced. This ion- 
isation front pushes a shock into the neutral gas, producing a 
sheared velocity field that is preferentially aligned with the x- 
axis (i. e., with the direction of the ionising photon field). This 
sheared velocity field eventually produces vortical motions that 
are perpendicular to both the x-axis and to the direction of the 
shear. This motion is seen in the xz-plane velocity field around 
the most massive clump in the t = 200 kyr frame shown in the 
bottom frame of Fig. 7. 

We have shown that the dense clumps that form as the re- 
sult of the photoevaporation of a dense, neutral structure in the 
ISM have angular momenta preferentially aligned in a direction 
perpendicular to the external ionising photon field. This result 
provides a natural explanation of the orientations observed in 



6 



V. Lora et al.: The angular momentum of condensations within elephant trunks 



the HH 555 (Bally & Reipurth 2003), HH 666 (Smith, Bally & 
Brooks 2004), and HH 333 (Rosado et al. 1999; Yusef-Zadeh 
et al. 2005) outflows, which emerge from elephant trunks in di- 
rections approximately perpendicular to the body of the trunks. 
Future observations of HH flows emerging from externally pho- 
toionised, neutral structures will show whether or not this kind 
of orientation is a general property of these outflows. 

We note again that we have simulated an ionisation front 
travelling into an initially steady, neutral medium with density 
perturbations. In this way, our simulations follow the dynam- 
ics produced by the propagating ionisation front and associated 
shock waves, which result in the production of clumps with an- 
gular momenta preferentially aligned perpendicular to the direc- 
tion of the ionising photon source. In the real ISM, a medium 
with density perturbations also has associated motions, and an 
initial vorticity field that will influence the angular momenta of 
clumps that might form (e. g., in the interaction with an ionisa- 
tion front). If the initial vorticity field is strong enough, it will 
probably hide the eff'ect of the vorticity generated by the shocks 
associated with the ionisation front, and the angular momentum 
alignment eff'ect described in this paper will not be present. 

An evaluation of whether or not the vorticity generated by 
the ionisation front will be hidden by the initial vorticity field of 
the cloud (present before the perturbations associated with the 
approaching ionisation front) can be completed on the basis of 
observations of the rotation of dense clumps in molecular clouds. 
For example, Ohashi et al. (1997) observed the kinematics of a 
number of NH3 cores in IRAS 04169-^2702 and computed their 
specific angular momenta. They find that cores with radii in the 
0.02 0.1 pc range have specific angular momenta L/M ~ 
(0.3 3)xlO^^ km^s~^ (clumps with larger radii having specific 
angular momenta up to an order of magnitude higher for a ~ 1 pc 
clump radius). 

In our simulation, the clumps with radii in the 0.02 ^ 0.1 pc 
range (corresponding to clump masses in the 0.1 100 M© 
range, see Sect. 3), have angular momenta L/M ~ (0.4 
20)x 10^^ km^s"^ (see Fig. 5). Therefore, our clumps have angu- 
lar momenta with values ranging from the lower L/M values of 
the cores observed by Ohashi et al. (1997), up to a factor of ~ 10 
times higher than the observed values. This result indicates that 
if the initial specific vorticity of the structure in the cloud were 
comparable to that of IRAS 04169-^2702, the passage of an ion- 
isation front would generate clumps of considerably higher spe- 
cific vorticity, and therefore the angular momentum alignment 
eff'ect described in this paper would indeed be present (at least 
for the more massive, higher angular momentum clumps). 

As a final point, we note that in the simulations presented 
in this paper we consider only the photoionisation of a neutral 
structure. In the case of the interaction of the radiation of an O 
star with a molecular cloud, it is unavoidable that the region out- 
side the ionisation front will be affected by the FUV radiation 
from the star, which at least partially photodissociates the ini- 
tially molecular material. Gorti & Hollenbach (2002) computed 
models of the photodissociation of dense clumps, and concluded 
that clumps with central column densities < 2 x 10^^ cm~^ (for 
an assumed cold-to-dissociated gas sound speed ratio of ~ 1/3) 
will be rapidly photodissociated, and disappear as local density 
enhancements. In our simulations, the clumps that are produced 
have central column densities of ~ (4, 9, 23, 47) x 10^^ cm"^ for 
clump masses of 0.01, 0.1, 1, and 100 M©, respectively (these 
central column densities are estimated by multiplying the clump 
radii given in Sect. 3 by the cutoff density of ~ 1.5 x 10^ cm~^). 
Therefore, in all cases the clumps have high enough column den- 
sities to avoid their dissipation by the incident FUV field. 



From the results of Gorti & Hollenbach (2002), we therefore 
conclude that the photodissociation caused by the FUV field will 
not destroy the clumps produced in our simulations. However, 
the early evolution of the flow (in which high density structures 
have not yet formed) might indeed be modified by the presence 
of a FUV field. It will therefore be interesting to carry out a 
future exploration of the formation of clumps within elephant 
trunks in the presence of both a photodissociating and a pho- 
toionising photon field. 
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